Interaction-round-a-face density- matrix 
renormalization-group method 

Wada Tatsuaki^'* and Nishino Tomotoshi^ 

^Department of Electrical and Electric Engineering, Faculty of Engineering, 
Ibaraki University, Hitachi 316-8511, Japan 

^Department of Physics, Faculty of Science, Kobe University, 657-8501, Japan 
Abstract 

We demonstrate the numerical superiority of the interaction-round-a-face (IRF) 
density-matrix renormahzation-group (DMRG) method apphed to SU (2) invariant 
quantum spin chains over the conventional DMRG. The ground state energy den- 
sities and the gap energies of both S = \ and S = 2 spin chains can be calculated 
using the IRF-DMRG without extensive computations. We have also studied the 
effect of tuning boundary interaction J^^^ at both ends of the chain from the IRF 
view point. It is clearly observed that the magnon distribution is uniform when the 
best Jgj^jj is chosen. 
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1 introduction 



Quantum spin chains are one of the most extensively studied systems by way 
of numerical calculations. In this article we consider isotropic spin-S* anti- 
ferromagnetic Heisenberg (AFH) spin chain, whose Hamiltonian is represented 
as 

Af-l 

H = JJ2S,- S,+i, (1) 

i=l 
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where J is a positive coupling constant and S,; denotes spin operator at i-th 
site. After Haldane conjectured the existence of finite excitation gap for integer 
spin chains and the absence for half-integer ones [3], many efforts have been 
performed to confirm this conjecture. We focus on numerical analysis of the 
AFH model in the following. 

A simple but concrete way to obtain the gap energy of the spin-S* AFH chain 
in the thermodynamic limit is to numerically calculate the excitation energies 
for finite size systems, and to perform the finite size scaling. Conventionally 
Lanczos method is applied to the Hamiltonian (Eq. (1)) expressed as a large 
sparse matrix in order to obtain lower lying states and their energies. The 
major shortcoming of this way is that the Lanczos diagonalization requires 
huge numerical resources. This is because the matrix dimension of the A^-site 
system is of the order of (25* + 1)^, that rapidly increases with both A^ and 
S. The largest system size N which is manageable with current computer is 
around 36 (for S — 1/2) even now. 

The Density-matrix renormalization-group (DMRG) invented by White [1,2] 
is a novel and powerful computational method, which overcomes the limitation 
of the manageable system size in the Lanczos diagonalization. It enables us to 
obtain physical observables such as ground state energy, correlation functions, 
etc., of large scale systems very accurately without performing extensive com- 
putations. It is not overwhelming to say that the DMRG is one of the standard 
tools to study low dimensional quantum systems. DMRG is a real-space renor- 
malization method, which divides the whole system into the left and the right 
blocks, and reduces the degree of freedom of each block using the eigenvalues 
of the reduced density matrix (DM) as probabilistic measures to choose the 
relevant block-spin states. The number of the block spin state kept under the 
RG transformation is conventionally represented by the letter m. It is impor- 
tant to choose sufficiently large m in order to keep numerical precision. The 
computational cost in DMRG is, roughly speaking, proportional to m^, which 
is not heavily dependent on size N. 

When applying DMRG to the spin-^ AFH chain specified by Eq. (1), the block 
spin states have total spin L up to nS, where n is a non-negative integer and 
depends on the number of spins in the block, we have to keep sufficient numbers 
of the block-spin state up to a certain total-spin L. This is problematic for large 
S chains, because the system has spin rotational symmetry; the eigenvalues 
of the DM are (2L + l)-fold degenerated, since a block spin state with spin L 
belongs to 2L-I-1 spin multiplet. Consequently the necessary block spin states 
(= m) increases rapidly with S. 

The problem of (2L -|- l)-fold degeneracy can be overcome by integrating out 
the SU (2) symmetry of the spin chain, and by representing the AFH Hamil- 
tonian using the total spin basis. As a result, the AFH system can be rep- 
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resented as a quantum limit of the interaction round a face (IRF) model [7]. 
Sierra and Nishino [4] have modified DMRG algorithm for the S = 1/2 and 
1 AFH models in the IRF representation (rep), where their numerical proce- 
dure is called 'IRF-DMRG'. Wada [6] has further developed IRF-DMRG for 
the cases S — 1 and 5" = 2. It is easily understood that the the computational 
cost of IRF-DMRG is much lower than that of the conventional DMRG, since 
the dimension of the associated Hamiltonian matrix in IRF rep is much smaller 
than that in conventional (= vertex) rep. 



2 IRF-DMRG applied to 6 = 1 and 2 AFH systems 



Now we demonstrate the numerical efficiency of IRF-DMRG when it is applied 
to AFH model for higher spin [S > 1) chains. In order to treat large S 
systems, Wada [6] has developed a way of obtaining the IRF expression for 
the Heisenberg interaction for 5" > 1. In the IRF formulation, the spin excited 
states can be targeted by putting an additional auxiliary (or ghost) spin to 
the spin chain; more details are shown in Ref. [6]. 

Before showing the calculated results we stress that all the numerical data 
are obtained using personal computer with the moderate power (SPECfp95 
~ 25). This is because the Hilbert space dimension in the IRF rep is much 
smaller than that in the conventional basis rep. For example, in the case 
of spin-2 spin chain, to keep miRp ~ 90 states in the IRF rep is equivalent to 
keep rug^ ~ 450 in the conventional rep. We keep 100 states at most, in the 
IRF rep. The ground state energy density Cq and the estimated gap energy A 
for 5 = 1 are eo = -1.401484039J and A = 0.4104J, and those for 5 = 2 are 
eo = -4.76124816 J and A = 0.088 J [6]. 

Having described the computational benefits of IRF formulation, now let us 

observe the boundary effect from the IRF view point. Since DMRG (including 
the IRF-DMRG) chiefly treats systems with open boundary, it is important 
to decrease the boundary effect and obtain the properties of bulk state. 

In order to estimate the spin excitation gap of AFH in the thermodynamic 
limit. White et al. [5] minimized the boundary effect by tuning the coupling 
constants Jcnd at both ends of the spin chain. They found that for S = 1 chain 
the best Jend is 0.5088 for the flrst excited state, and is 0.7 for the ground 
state. Let us review how fast the boundary effect decreases with the system 
size when Jend = 0.5088. Figure 1 shows the gap energies A as a function of 
the system size N . After rapid increase with iV, the gap energy A(iV) becomes 
almost independent of N when it is larger than 20. It should be noted that it 
is possible to estimate A = 0.4015J from the data < 20. In contrast, when 
Jend deviates from the best value, A(A^) shows very slow increases with N, 
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this makes the finite size scahng very difficult. 

What happens when Jgnd = 0.5088 for the first excited state? It is believed 
that the magnon (or the domain wall excitation) has zero momentum (k — n 
in their definition) and it distributes quite uniformly [5]. The situation is 
quite similar to the case of a particle in tight binding model with shallow 
negative potential at the system boundary. The IRF rep detects the magnon 
distribution more clearly, since it is possible to obtain the partial spin average 

{Su)^d:s,). (2) 

The derivative of {3^) with respect to the lattice coordinate i is proportional 
to the density of the quasi particle (= magnon). To speak precisely, it is 
possible to obtain {S^) from the numerical data of the conventional DMRG, 
but careful treatment of numerical data is required; the IRF-DMRG provides 
(Sij) more directly. Figure 2 shows (S^^) for three different Jend s. As is seen 
the partial spin average is quite linear in i at the best boundary condition 
(BC) of Jcnd = 0.5088. The partial spin average shall be a good reference 
data when searching the similar BC for the higher spin chains. More detailed 
analysis will be pubhshed elsewhere. 

In conclusion we have demonstrated the numerical advantage of IRF for- 
mulation through the study of the AFH model, particularly for higher spin 
chains. The ground state energy densities Cq and gap energies A are obtained 
with moderate computational cost. The results are Cq — — 1.401484039J, 
A = 0.4104J for S ^ 1, and Cq = -4.76124816J, A = 0.088J for S = 2. 
The estimated spin gap is much more accurate than those obtained so far, 
under the condition m < 100. In addition, we have discussed how the min- 
imization of the boundary effect is observed from the IRF view point. We 
have seen that at the best boundary condition at which Jend minimizes the 
boundary effect, the partial spin average {3^^^} = Sj) of the first excited 

state is almost linear in i. To find out such smooth BCs for higher spin cases 
is our future subject. 

T. N thank to G. Sierra for valuable discussions. The present work is supported 
by Grant-in- Aid for Scientific Research from Ministry of Education, Science, 
Sports and Culture (No. 11640376). Most of computations were performed on 
a 600 MHz Alpha PC164 workstation with a variant of Linux Alpha OS. 
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Fig. 1. The gap energies A as a function of the chain length N for spin-1 AFH spin 
chain. The number of states m being kept during IRF- DMRG iterations are only 
32. The Jend denotes the coupling constants at both ends of the chain. All coupling 
constants J except Jend set to unity. When Jend sets to 0.5088, A{N) is almost 
independent of N. 



1.2 




Q 5 I I I I I I I I I I I 

10 20 30 40 50 60 70 80 90 100 



site index i 



Fig. 2. The partial spin average {Su) distribution for the first excited state of the 
spin-1 AFH chain. When Jend = 0.5088, at which the boundary effect is minimized, 
{Sii) becomes linear in site index i. 
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